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EXECUTIVE  SUMMARY 


A  fuel  generation  model  for  char  forming  polymers  in  fires  is  developed  from  mechanistic 
pyrolysis  kinetics.  Under  conditions  of  flaming  combustion  the  coupled  rate  equations  for  thermal 
degradation  products  and  reactants  reduce  to  a  single  rate  law  for  the  residual  mass.  Exact  results 
are  obtained  from  the  mass  loss  history  which  include  an  equilibrium  char  yield  whose  value 
depends  only  on  the  relative  rates  of  gas  and  char  formation  at  a  particular  temperature.  Reaction 
rate  constants  for  thermolysis  of  chemical  bonds,  gas  production,  and  char  formation  are 
determined  from  parametric  fits  of  the  mechanistic  charring  model  to  thermogravimetric  data. 
Predictions  of  the  nonisothermal  mass  loss  during  constant  heating  rate  experiments  are  in 
agreement  with  experimental  data  over  the  expected  range  of  validity. 
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INTRODUCTION 


Full-scale  fire  tests  have  demonstrated  that  the  heat  release  rate  of  a  burning  material  is  the  single 
most  important  parameter  determining  its  fire  hazard  in  an  enclosed  environment  such  as  a  room 
[1],  submarine  [2],  rail  car  [3],  or  passenger  aircraft  cabin  [4],  Bench-scale  fire  tests  have  shown 
that  polymers  which  form  a  carbonaceous  char  during  burning  have  lower  ignitability  [5]  and 
lower  heat  release  rate  than  nonchar  forming  polymers  [6].  Char  formation  in  a  fire  limits  the 
amount  of  volatile  fuel  that  can  be  produced  by  the  burning  polymer,  provides  a  thermally 
insulating  layer  at  the  surface  to  reduce  heat  transmission  into  the  material,  and  acts  as  a  diffusion 
barrier  to  combustible  gases.  Despite  the  empirical  correlation  between  char  formation  and  reduced 
flammability  of  synthetic  polymers  [5-7],  pyrolysis  models  which  include  char  formation  have 
focused  almost  entirely  on  the  naturally  occurring  polymer  cellulose  [8-19],  The  present  work 
provides  a  generalized  pyrolysis  model  which  accounts  for  thermally  induced  gasification  and  char 
formation  of  synthetic  polymers  and  which  can  be  evaluated  using  standard  laboratory 
thermogravimetric  equipment.  Analytic  solutions  having  a  minimum  number  of  adjustable 
parameters  were  pursued  in  order  to  provide  insight  into  the  solid  state  thermochemical  processes 
governing  polymer  combustion  and  to  help  guide  the  search  for  fire  resistant  materials  for  next 
generation  aircraft  interiors. 


BACKGROUND 

The  steps  involved  in  the  flaming  combustion  of  solid  polymers  are  described  schematically  in 
figure  1,  after  van  Krevelen  [20],  Flaming  combustion  requires  three  coupled  processes:  heating 
of  the  polymer,  thermal  decomposition/pyrolysis,  and  ignition  of  the  gaseous  decomposition 
products  in  air.  An  ignition  source  or  thermal  feedback  of  radiant  energy  from  the  flame  supplies 
heat  to  the  polymer  surface  which  causes  thermolytic  cleavage  of  primary  chemical  bonds  in  the 
polymer  molecules.  Pyrolysis  gases  mix  and  react  with  air  in  the  combustion  zone  above  the 
surface  releasing  heat  and  producing  carbon  dioxide,  water,  and  incomplete  combustion  products 
such  as  carbon  monoxide  and  soot. 
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FIGURE  1.  SCHEMATIC  OF  THE  COMBUSTION  PROCESS  OF  SOLID  POLYMERS 

AFTER  [20] 

One-dimensional  burning  models  for  polymers  have  been  developed  with  coupled  heat  and  mass 
transfer  and  kinetically  controlled  thermal  decomposition  in  the  pyrolysis  zone  of  the  burning 
material  for  both  chairing  [16-17]  and  noncharring  polymers  [21],  The  basic  thermal  degradation 
mechanism  leading  to  volatile  fuel  generation  in  char  forming  polymers  has  been  described  as  a 
generalized  chemical  bond  scission  process  consisting  of  primary  and  secondary  decomposition 
events  [19,  22-24],  The  primary  decomposition  step  is  main-,  end-,  or  side-chain  scission  of  the 
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polymer  to  form  free  radical  intermediates.  Subsequent  hydrogen  transfer  and/or  recombination  of 
the  intermediates  lead  to  primary  volatiles  (gas  and  tar)  and  char.  The  primary  char  further 
decomposes  by  dehydrogenation  to  form  the  secondary  gas  and  a  thermally  stable  secondary  char. 

Broido  and  Nelson  [9]  may  have  been  the  first  to  propose  a  mechanistic  parallel  reaction  model 
wherein  cellulose  thermally  decomposes  directly  to  tar  and  char  through  competing  reactions.  The 
two-step  parallel  reaction  model  for  cellulose  predicted  a  variable  but  finite  residual  mass  after 
pyrolysis  which  depended  on  time  and  temperature  (i.e.,  heating  rate)  through  the  magnitude  of  the 
rate  constants  for  the  competing  processes.  Lewellen  et  al.  [13]  extended  the  cellulose  parallel 
reaction  scheme  to  include  the  formation  and  evaporation  of  an  active  intermediate  (levoglucosan) 
which  subsequently  degrades  to  tar  and  char  through  competing  reactions.  The  introduction  of  a 
reactive  intermediate  by  Lewellen  qualitatively  accounted  for  the  absence  of  char  in  rapidly  heated 
(102— 105  K/s)  samples  where  the  residence  time  of  the  intermediate  in  the  condensed  phase  was 
insufficient  to  allow  char  formation,  i.e.,  the  rate  of  char  formation  was  slow  compared  to  the  rate 
of  evaporation  of  the  reactive  intermediate  at  high  temperatures.  Suuberg  et  al.  [15]  adopted 
Lewellyn’s  model  of  cellulose  decomposition  but  suggested  that  the  first  step  in  the  pyrolysis- 
generation  of  the  reactive  intermediate  from  the  parent  cellulose—was  a  reversible  process,  as 
suggested  for  levoglucosan  formation  [25],  The  reversibility  of  the  first  step  implied  a  low 
concentration  of  the  reactive  intermediate  during  the  course  of  the  pyrolysis  which  accounted  tor 
the  low  measured  vapor  pressures.  Moreover,  equilibrium  between  the  cellulose  and  tar 
intermediate  would  not  be  achieved  at  low  temperatures  because  the  intermediate  would  be  rapidly 
consumed  by  gasification  and  chairing  [15].  Reversible  dissociation  of  the  polymer  to  a  reactive 
intermediate  is  analogous  to  the  ceiling  temperature  in  polymerization  reactions  [26]. 

Suuberg’s  model  suggests  that  formation  of  the  active  tar  intermediate  is  the  rate  limiting  step  at 
temperatures  below  the  tar  vaporization  temperature  =  530-570  K  [15,  27],  while  mass  transport  of 
the  tars  is  the  rate  limiting  step  at  high  heating  rates  and  temperatures.  A  schematic  representation 
of  Suuberg’s  model  is 


cellulose  _^T  tar  intermediate 


evaporated  tar  +  some  gas 
char  +  light  gases 


Decomposition  schemes  which  account  for  some  or  all  of  these  pyrolysis  products  have  been 
proposed  wherein  the  decomposition  steps  occur  sequentially  (series),  simultaneously  (parallel),  or 
in  some  combination  of  series/parallel  steps  [9-19].  Three  of  these  mechanistic  pyrolytic  reaction 
schemes  for  cellulose  have  been  reviewed  recently  [18]  including  the  single-step  first  order  model 
[17],  an  uncoupled  three-step  parallel  model  [8],  and  a  coupled  three-step  series-parallel  model  [9, 
19].  The  single-step  and  three-step  uncoupled  models  have  a  fixed  char  yield  as  an  adjustable 
parameter  while  the  three-step  coupled  model  has  a  variable  char  yield.  All  of  the  models  predict 
rate  dependent  peak  decomposition  temperatures.  Variable  (n-th)  order  decomposition  kinetics 
have  been  fit  to  mass  loss  data  for  char  forming  polymers  with  reasonable  success  using  reaction 
order  [12,  28]  and  empirical  weighting  factors  [29]  as  adjustable  parameters.  However,  little 
insight  is  gained  into  the  reaction  pathways  from  these  curve-fitting  exercises.  An  analytic  solution 
for  the  isothermal  residual  mass  of  the  polymer  in  the  single-step  and  three-step  uncoupled  and  the 
n-th  order  empirical  models  is  possible  because  the  reactions  are  first  order  in  the  polymer  mass. 
The  coupled  model  has  been  solved  analytically  for  the  residual  isothermal  mass  but  requires  six 
kinetic  parameters  and  a  maximum  char  yield  for  evaluation.  Numerical  solutions  of  the  three-step 
delinked  [16]  and  parallel  models  have  been  published  for  constant  and  variable  heating  rates  [18]. 


2 


POLYMER  THERMAL  DECOMPOSITION 


The  condensed  phase  thermal  decomposition  pathway  of  polymers  is  assumed  to  be  a  series- 
parallel  reaction  model  [9,  19]  which  includes  a  reversible  first  step  [15]  and  secondary  char 
formation  [20]  as  shown  in  figure  2. 
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FIGURE  2.  SCHEMATIC  REPRESENTATION  OF  PRIMARY  AND  SECONDARY 
DECOMPOSITION  PROCESSES  IN  POLYMER  PYROLYSIS 

Equations  1-5  describe  the  anaerobic  reactions,  and  equations  6-8  extend  the  model  to  include 
condensed  phase  oxidation  reactions  of  the  reactive  intermediate,  primary,  and  secondary  char. 
Bulk  and  surface  oxidation  reactions  become  important  for  long-term  thermoxidative  stability  and 
smoldering  combustion,  respectively.  In  reactions  1-8,  P  is  the  polymer,  I*  is  the  reactive 
intermediate,  G1  and  Cl  are  the  primary  volatiles  and  char,  and  G2  andC2  are  the  secondary  gas 
and  char.  Solid-state  bulk  and  surface  oxidation  reaction  products  for  the  intermediate,  primary 
char,  and  secondary  char  are  denoted,  UO2,  CIO2,  and  C2O2 ,  respectively. 
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In  equations  1-8,  kp  and  k_p  are  the  rate  constants  for  the  forward  and  reverse  polymer  dissociation 
steps,  respectively;  kgi  and  kd  are  rate  constants  for  primary  gas/tar  and  char  formation;  kg2  and 
kC2  are  the  rate  constants  for  secondary  gas  and  char  formation;  and  kpo,  kdo,  and  k^o  are  the 
rate  constants  for  the  oxidation  reactions  of  the  intermediate  species  and  primary  and  secondary 
chars  in  the  solid  state.  Although  equations  1-8  represent  a  highly  simplified  set  of  reactions  and 
products,  the  resulting  system  of  rate  equations  is  highly  coupled  and  an  analytic  solution  in  terms 
of  the  overall  sample  mass  is  impossible. 

MASS  LOSS  KINETICS 

The  objective  of  this  work  is  to  derive  analytic  solutions  for  the  mass  loss  history  of  a  burning 
polymer  from  a  mechanistically  based  kinetic  model  of  polymer  pyrolysis  which  accounts  for 
products  that  are  important  in  the  burning  process,  i.e.,  combustible  gases  and  char.  Of  particular 
interest  are  simple  solutions  for  the  mass  loss  history  of  a  polymer  which  provide  insight  into  the 
mechanisms  of  decomposition,  allow  estimation  of  the  fuel  generation  rate  under  isothermal  and 
nonisothermal  heating,  and  can  be  verified  using  standard  laboratory  thermogravimetric 
techniques.  To  develop  the  model,  assumptions  are  made  about  the  process  of  polymer  thermal 
degradation  which  are  relevant  to  flaming  combustion  and  which  will  reduce  equations  1-8  to  a 
tractable  set.  The  following  are  proposed: 

Assumption  A:  The  breaking  of  primary  chemical  bonds  in  the  polymer  (thermolysis)  [30]  is  the 
rate  limiting  step  at  heating  rates  and  surface  temperatures  observed  in  fires  [31],  i.e.,  P  ->  I* 
is  slow  compared  to  the  I*  — >  G1  and  I*  — >  Cl  reactions  which  form  volatiles  and  char. 

Assumption  B:  The  reactive  intermediate  I*  generated  in  the  polymer  dissociation  step  is  in 
dynamic  equilibrium  with  the  parent  polymer  but  is  consumed  in  the  process  of  gas  and  char 
formation  such  that  its  concentration  never  becomes  appreciable  and  decreases  slowly  over  time 
as  the  polymer  is  consumed.  Consequently,  the  rate  of  change  of  I*  with  time  is  insignificant 
compared  to  the  rate  of  polymer  consumption,  gas  production,  and  char  formation,  so  that  for 
computational  purposes,  dl*/dt  =  0.  This  is  the  stationary-state  hypothesis. 

Assumption  C:  Thermal  decomposition  of  primaiy  char  to  secondary  char  and  gas,  i.e..  Cl  ->  C2 
and  Cl  — >  G2,  is  slow  compared  to  the  formation  of  the  primary  char  I*  — >  Cl  at  typical 
flaming  surface  temperatures  of  350-700°C.  Consequently,  only  the  primary  char  is 
considered  in  formulating  the  reaction  set  for  the  mass  loss  model. 

Assumption  D:  The  thermoxidative  environment  in  the  pyrolysis  zone  of  a  burning  solid  polymer 
is  anaerobic.  Dissolved  molecular  oxygen  and  oxygen  diffusion  into  the  pyrolysis  zone  of  the 
solid  are  considered  negligible  with  respect  to  their  effects  on  gas  and  char  formation  so  that 
equations  6-8  can  be  neglected  in  the  kinetic  model  of  polymer  combustion.  This  assumption 
does  not  preclude  the  possibility  of  surface  mass  loss  due  to  thermoxidative  reactions  at  the 
polymer-air  interface,  e.g.,  thermogravimetric  (TGA)  experiments  conducted  in  air  or 
smoldering  fires. 
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In  combination  with  the  generalized  combustion  and  pyrolysis  schemes  of  figures  1  and  2, 
assumptions  A-D  lead  to  a  simplified  mass  loss  model  for  polymer  combustion  which  is  shown  in 
figure  3.  This  simplified  scheme  reduces  gas  and  char  formation  to  a  single  step  involving  parallel 
reactions  of  the  active  intermediate,  I*. 
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FIGURE  3.  KINETIC  MODEL  FOR  POLYMER  FLAMING  COMBUSTION 
The  reduced  reaction  set  from  figure  3  becomes: 
k„ 


I* 


-p 


I*  - G 


kc 


I* — C 

and  the  system  of  rate  equations  for  the  species  at  time  t  is 

f  =  -  V  +k-p!* 


dl* 

dt 


kpP-(k_p  +  kg+kc)I* 


dG  _  u  t* 

dt  ~  s 

dC  —  lc  U 
dt  ~  Kcl 

According  to  the  stationary- state  hypothesis  (Assumption  B),  dl*/dt  =  0,  so  that  equation  12 
provides  the  useful  result 


(1) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 
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p  = 


P  =  k*P 


[k_p  +  kg  +  kcJ 

Substituting  P  =  k*P  from  equation  15  into  equations  11, 13,  and  14 


(15) 


dP 

dt 


=  -[kp-k*k-P]P 


lit  =  k*k*  p 
f  =  p 

With  P  «  P,  G,  and  C  the  total  mass  balance  in  terms  of  the  initial  mass,  m0,  is: 

m0  =  P  +  G  +  C  +  P  =  P  +  G  +  C 

Since  dm0/dt  =  0,  equations  16  and  19  lead  to 

dP  _  dC  dG  _  _ [u  _ Ic^lc  IP 

dt  ~  dt  dt  ~  l-kP  k-PJ P 


(16) 

(17) 

(18) 

(19) 

(20) 


The  sensible  mass,  m,  of  the  sample  as  measured,  for  example,  in  a  TGA  experiment  or  fire 
calorimeter  test  is 


m  =  P  +  C  +  P  =  P  +  C 
From  equations  17  and  20 


dm  _  c/P  ,  dC  _  _  dG  _  _  u*),  p 
dt  ~  dt  +  dt  ~  dt  ~  s 


(21) 


(22) 


Equation  16  can  be  solved  immediately  for  P  for  the  initial  condition,  P  =  P0  =  m0  @  t  =  0,  with 
the  result 


P  =  m0exp(-[kp-k*k_p]t) 

Substituting  this  result  for  P  into  equation  22  and  separating  variables 

r  m  r  t 

dm’  =  -  k+kgm0  exp  (-  [kp  -  k*k_p]  t)  dt 
Jo 


(23) 


(24) 
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Since  kg  and  kc  »  kp  and  k_p  by  Assumption  A, 


the  solution  of  equation  24  is 


kp-k*k_p  =  kp, 


kg  +  kc 


k  +k  +kc 


m(t) 

m0 


1- 


kg 


kg  +  k£ 


{1  -  exp  (-kp  t)} 


(25) 


Equation  25  shows  that  as  t  — >  °°  the  residual  mass,  m/m0,  approaches  an  equilibrium  value  at 
constant  temperature  given  by 


m(<») 


m, 


Lkg+kc 


=  Y„ 


(26A) 


Assuming  Arrhenius  forms  for  kg  and  kc, 

A 


YC(T) 


l-i 


1  +  ^exp[-(Eg-Ec)/RT] 


(26B) 


with  Ec  and  Eg  and  Ac  and  Ag  being  the  activation  energies  and  frequency  factors  for  char  and  gas 
formation,  respectively.  Equation  26  predicts  a  finite  residual  mass  at  infinite  time  if  kc  >  0  and 
zero  residual  mass  if  kc  =  0.  The  physical  significance  of  a  temperature  dependent,  equilibrium, 
residual  mass  as  the  ratio  of  rate  constants  for  volatile  and  char  formation  is  consistent  with  the  use 
of  group  contributions  for  the  char  forming  tendency  of  polymers  at  a  particular  temperature  [20]. 


The  crossover  temperature,  Tcr,  is  defined  [18]  as  the  temperature  at  which  the  rates  of  gas  and 
char  formation  are  equal,  i.e.,  when  kg  =  kg.  From  equation  26B 


T, 


er 


(Eg~Ec) 
R  ln[Ag/Ac] 


(27) 


It  follows  from  equation  26A  that  the  crossover  condition,  kg  =  kc,  corresponds  to  an  equilibrium 
char  yield  of  50%,  i.e.,  Yc  =  0.50. 


Substituting  equation  26A  into  equation  25  yields  the  final  result  for  the  isothermal  mass  loss 
history  in  terms  of  the  rate  constant  kp  for  thermolysis  of  primary  chemical  bonds  in  the  polymer. 


W;  =  Yc  +  (l-Yc)exp(-kp.) 


(28) 


If  Yc  =  ji  =  constant,  equation  28  is  the  solution  for  the  isothermal  mass  loss  history  of  a  filled 
polymer  with  a  nonvolatile  mass  fraction  |i  satisfying  the  rate  law 

=  ~  kp  (  m  -  (im0 )  (29) 

although  equation  29  was  not  assumed  a  priori  in  the  present  derivation. 
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While  estimation  of  the  material  parameters  governing  polymer  pyrolysis  is  best  accomplished 
under  isothermal  conditions  as  above,  many  processes  of  interest  such  as  the  production  of 
combustible  gases  in  a  fire  occur  at  finite  heating  rates.  Consequently,  it  is  important  to  be  able  to 
predict  mass  loss  under  dynamic  heating  conditions.  For  a  constant  heating  rate,  {3  =  dT/dt,  the 
independent  variable  in  equation  28  can  be  transformed  from  time  to  the  dimensionless  variable 


x 


kp(T)  t(T)  =  - 


A  (T  -  T0 ) 

P 


exp 


(30) 


so  that  equation  28  can  be  written  in  the  elementary  form 


=  (l-Yc)fXmaVdx  (31) 

O  Jo 

The  exact  solution  of  equation  31  for  the  fractional  mass  as  a  function  of  temperature  at  constant 
heating  rate  is 


-^(T-T0)e-1VRT]  (32) 

Equation  32  is  within  0.001%  of  the  value  obtained  by  double  precision  numerical  integration  of 
equation  31  over  the  temperature  interval  T0  =  273  K,  T  =  1273  K  at  1  K  increments  after 
appropriate  transformation  of  the  independent  variable  from  x  to  T. 

The  fractional  mass  loss  rate  during  a  linear  temperature  program  is  obtained  by  differentiating 
equation  32  with  respect  to  time.  No  simple  solution  is  possible  when  Yc  =  YC(T)  so  we  make  the 
drastic  approximation  that  the  equilibrium  residual  mass  (char  yield)  is  constant  and  independent  of 
temperature,  i.e.,  Yc  -  p  and  obtain 


mu; 

m„ 


=  Yc  +  (1-YJexp 


-1  dm  „ 
dt 


(1  -\x)kp(T) 


,  .  (T-T„)E„\ 

1 + ~w~  n 


(T-TJkJT) 


P 


(33) 


Noting  that  the  maximum  rate  of  change  of  m/m0  versus  time/temperature  occurs  when  x  =  -1,  we 
set  the  exponent  in  equation  33  to  -1  and  obtain  the  useful  approximation  for  the  maximum  mass 
loss  rate  of  a  polymer  with  inert  fraction  p  at  constant  heating  rate 


-1  dm 
dt 


(1-M) 


PEa 

eRT,:,. 


(34) 


where  Tmax  is  the  temperature  at  peak  mass  loss  rate  obtained  from  the  root  of  [32] 


Ea  1 

2 

1 

Ea  1 

4.  In 

RP] 

.RTmux. 

I 

,RTmax. 

i  111 

[AEaJ 

(35) 
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EXPERIMENTAL 


Anaerobic  pyrolysis  of  3-5  mg  samples  of  low-density  polyethylene  (Mw  -  35,000  g  mol-1, 
Mw/Mn  =  4.5,  and  p  =  906  kg  m~3,  Aldrich  Chemical),  polymethylmethacrylate  (PMMA,  Aldrich 
Chemical,  Mw  =  15,000  g-moH)  and  crosslinked  phenolic  triazine  thermoset  resin  (PRIMASET 
PT  30,  Lonza  Chemical,  or  AROCY  XU  371,  Ciba  Specialty  Chemicals)  were  obtained  under 
nitrogen  (99.998%  N2,  Matheson)  flowing  at  0.10  L  min^1  in  a  therm ogravimetric  analyzer 
(Perkin  Elmer  TGA-7)  after  a  15-minute  purge  to  remove  residual  air.  Anaerobic,  isothermal  mass 
loss  histories  were  recorded  at  temperatures  ranging  from  350-450°C  for  the  phenolic  triazine  and 
nonisothermal  scans  were  performed  at  selected  constant  heating  rates  ranging  from  1  to  100 
K/min  for  all  materials.  The  polymethylmethacrylate  and  low-density  polyethylene  were  used  as 
received.  The  high-purity  trifunctional  cyanate  ester  novolac  monomer  was  used  as  received 
without  further  purification  and  thermally  cured  without  catalyst  for  four  hours  at  250°C  to  >  95% 
conversion  as  determined  by  infrared  spectroscopy  [33],  Curing  the  cyanate  ester  novolac 
monomer  produces  a  tightly  crosslinked  thermoset  network  of  oxygen-linked  triazine  rings 
(cyanurates)  having  the  repeat  unit  atomic  composition  QH5NO  and  the  chemical  structure 


Experimental  data  for  the  isothermal  mass  loss  histories  of  the  cured  phenolic  triazine  resin  are 
shown  as  solid  symbols  in  figure  4.  Solid  lines  are  the  best  fit  of  equation  28  to  the  experimental 
data  using  kp  and  Yc  as  adjustable  parameters.  Table  1  lists  values  for  kp  and  Yc  obtained  from  the 
fit  of  equation  28  to  experimental  data  in  figure  4.  Figure  5  is  an  Arrhenius  plot:  ln(kp)  =  ln(A)  - 
Ea/RT  of  the  T  and  kp  values  in  table  1.  The  slope  in  figure  5  is  Ea/R  from  which  the  activation 
energy  is  Ea  =  178  kJ/mol.  The  intercept  is  the  natural  logarithm  of  the  frequency  factor  for 
pyrolysis  of  the  phenolic  triazine,  A  =  109  sec1. 
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FIGURE  4.  MEASURED  AND  CALCULATED  ISOTHERMAL  MASS  LOSS  HISTORIES  AT 
INDICATED  TEMPERATURES  FOR  PHENOLIC  TRIAZINE  RESIN.  (Solid  lines  are  fit  of 
equation  28  using  values  for  kp  and  Yc  in  table  1.) 

It  is  important  to  recognize  that  both  the  activation  energy  and  the  frequency  factor  for  pyrolysis 
obtained  from  the  mass  loss  model  for  char  forming  polymers  (28)  are  significantly  lower  than  the 
values  for  these  parameters  obtained  using  generalized  models  for  analysis  of  thermogiavimetric 
data.  For  example,  the  Arrhenius  parameters  for  kp  of  the  phenolic  triazine  resin  obtained  in  this 
work  (Ea  =  178  kJ/mol,  A  =  1()9  sec1)  compare  to  Ea  =  280  -  326  kJ/mol  and  A  =  1018  -  10:0 
sec1  obtained  from  the  same  data  using  n-th  order  models  [32,  33]. 

TABLE  1.  BEST  FIT  kP  AND  Yc  VALUES  FOR  PHENOLIC  TRIAZINE  THERMOSET  RESIN 


T(°C) 

kn  (sec-1) 

Yc 

350 

1.4  x  1()-5 

0.92 

400 

8.5  x  1()-5 

0.90 

420 

2.7  x  1(H 

0.88 

450 

1.2  x  1(L3 

0.84 

The  relative  rate  constants  for  volatile  and  char  formation  are  obtained  by  plotting  equation  26B  in 
the  form 


In 


1-YC1 

=  In 

rA8l 

Eg  ~  Ec" 

Ye  J 

[AcJ 

R  J 

1 

T 


(36) 
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FIGURE  5.  PLOT  OF  THE  NATURAL  LOGARITHM  OF  kP  VERSUS  RECIPROCAL 
TEMPERATURE  FOR  PHENOLIC  TRIAZINE  RESIN 


Thus,  a  plot  of  ln[(l-Yc)/Yc]  versus  1/T  has  a  slope  proportional  to  the  difference  in  activation 
energies  for  volatile  and  char  formation  and  an  intercept  which  is  the  natural  logarithm  of  the  ratio 
of  the  frequency  factors.  Figure  6  shows  a  plot  of  equation  36  using  the  isothermal  char  yields  in 
table  1.  The  slope  gives  Eg  -  Ec  =  +30  kJ/mol  and  the  intercept  Ag/Ac  =  17.  This  result  indicates 
that  even  though  the  molecular  collision  frequency  for  gas  production  is  seventeen  times  greater 
than  for  char  formation,  char  formation  has  a  lower  activation  energy  by  30  kJ/mol  which  is 
perhaps  due  to  resonance  stabilization  of  the  polycyclic  aromatic  char  [34].  The  Arrhenius 
parameters  for  gas  and  char  formation  substituted  into  equation  26B  allow  calculation  of  the 
equilibrium  char  yield  as  a  function  of  temperature.  The  crossover  temperature  for  the  phenolic 
triazine  resin  from  equation  27  is 


(30  kJ/mol) 
R ln[ 17] 


1273  K 


For  the  phenolic  triazine  resin  the  rates  of  primary  char  and  volatile  formation  are  equal  at  1000°C 
and  Yc=  0.50.  In  comparison  the  crossover  temperature  for  cellulose  is  about  250°C  from  char 
yield  data  and  the  criteria  Yc  =  0.50  [9]. 
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FIGURE  6.  PLOT  OF  THE  NATURAL  LOGARITHM  OF  (1-YC)/YC  VERSUS  RECIPROCAL 
TEMPERATURE  FOR  PHENOLIC  TRIAZINE  RESIN 

The  predicted  mass  loss  of  the  phenolic  triazine  resin  at  constant  heating  rates  of  (3  =  1,  5,  and  20 
K/min  calculated  from  equation  32  is  compared  to  thermogravimetric  data  in  figure  7.  The 
parametrically  determined  A  =  109  sec*1,  Ea  =  178  kJ/mol,  and  the  temperature  dependent  char 
yield  YC(T)  from  equation  26B  were  used  in  the  calculations.  Solid  symbols  are  the 
thermogravimetric  data  at  the  indicated  heating  rates.  Solid  lines  are  the  predictions  of  equation  32 
for  those  same  heating  rates,  and  the  dotted  line  is  the  equilibrium  char  yield  corresponding  to 

P  =  0. 


Agreement  between  the  measured  and  calculated  residual  mass  of  the  phenolic  triazine  resin  in  the 
scanning  TGA  experiments  is  generally  good  for  temperatures  below  about  700°C  as  shown  in 
figure  7.  The  difference  between  the  model  and  the  experimental  data  above  700°C  at  these  heating 
rates  is  a  result  of  neglecting  secondary  charring  and  gasification  in  the  simplified  mass  loss  model 
which  considered  only  primary  char  and  gas  formation  (Assumption  C). 

Equation  32  was  used  to  predict  mass  loss  curves  for  the  phenolic  triazine  resin  at  the  high  heating 
rates  (3  =  1  to  100  K/s  characteristic  of  a  burning  polymer  in  a  fire  [31].  Figure  8  shows  these 
calculations  for  five  decades  of  heating  rate  ranging  from  (3  =  0.01  to  100  K/s.  Figure  8  shows 
that  both  the  temperature  at  maximum  mass  loss  rate,  Tmax,  and  the  volatile  mass  fraction,  m/m0  - 
Yc  (Tmax),  increase  significantly  with  heating  rate  for  a  chairing  polymer  as  the  material  becomes 
superheated. 
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Temperature  (°C) 


FIGURE  7.  COMPARISON  OF  MEASURED  AND  CALCULATED  MASS  LOSS  OF 
PHENOLIC  TRIAZINE  RESIN  IN  THERMOGRAVIMETRIC  ANALYZER  AT  CONSTANT 
HEATING  RATES  OF  (3=1,5,  AND  20  K/min.  (Solid  lines  are  fit  of  equation  32  to 
experimental  data  for  indicated  heating  rates.  Dashed  line  is  calculated  equilibrium  char  yield 

(p  =  0)  versus  temperature.) 


FIGURE  8.  MASS  LOSS  CALCULATIONS  FOR  PHENOLIC  TRIAZINE  RESIN  AT  FIVE 

DECADES  OF  HEATING  RATE 
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Table  2  compares  measured  and  calculated  peak  fractional  mass  loss  rates  at  various  heating  rates 
for  the  noncharring  polymers  PMMA  and  polyethylene  and  the  char-forming  phenolic  triazine 
resin.  Calculated  values  were  obtained  using  equation  34  with  the  indicated  heating  rates.  For  the 
phenolic  triazine  with  Ea  =  178  kJ/mol,  a  value  of  |i  =  0.75  =  Yc  (650°C)  gave  the  best  fit  of  the 
peak  mass  loss  rate  data.  For  the  polymethylmethacrylate  of  this  study,  Ea  =  160  kJ/mol  [32,  35] 
and  ft  =  0,  while  for  low-density  polyethylene,  Ea  =  200  kJ/mol  [35,  36]  and  ji  =  0.  The  average 
error  in  the  calculated  peak  mass  loss  rates  is  about  ±  30%  of  experimental  values  which  is  on  the 
order  of  the  uncertainty  in  the  activation  energies  used  in  the  calculation.  Thus,  equation  34 
qualitatively  describes  the  peak  mass  loss  rate  of  charring  and  noncharring  polymers  over  a  wide 
range  of  constant  heating  rates.  This  result  could  explain  the  lower  heat  release  rate  of  polymers 
with  high  decomposition  temperature  and  char  yield  in  terms  of  the  mass  loss  (fuel  generation)  rate 
in  the  pyrolysis  zone  during  flaming  combustion. 

TABLE  2.  MEASURED  AND  CALCULATED  PEAK  MASS  LOSS  RATES  FOR 
POLYMETHYLMETHACRYLATE,  POLYETHYLENE,  AND  PHENOLIC  TRIAZINE 

POLYMERS 


Heating 

Rate 

Tmax 

(measured) 

Peak  Mass  Loss 
Rate  (measured) 

Peak  Mass  Loss 
Rate  (calculated) 

K/min 

Kelvin 

mg/g-s 

mg/gs 

PHENOL 

IC  TRIAZINE 

1 

693 

0.03 

0.07 

3 

703 

0.14 

0.20 

5 

724 

0.29 

0.31 

10 

743 

0.63 

0.59 

20 

745 

1.28 

1.18 

POLYMETHYLMETHACRYLATE 

3 

639 

1.38 

1.03 

10 

653 

3.81 

3.32 

30 

682 

10.03 

9.65 

100 

698 

28.63 

31.11 

200 

719 

47.65 

61.03 

POLYETHYLENE 

3 

738 

1.89 

0.81 

10 

757 

5.55 

2.63 

30 

778 

14.27 

7.31 

100 

790 

26.38 

23.68 

200 

808 

63.50 

45.14 

CONCLUSIONS 

The  objective  of  this  work  was  to  obtain  analytic  solutions  for  the  mass  loss  and  mass  loss  rate  of 
char  forming  polymers  from  a  mechanistic  pyrolysis  model  to  provide  insight  into  the  fuel 
generation  process  during  flaming  combustion.  Analytic  solutions  for  the  residual  mass  having  a 
minimum  number  of  adjustable  parameters  were  desired  so  that  the  model  could  be  easily  evaluated 
for  different  materials  on  conventional  therm ogravimetric  equipment.  A  mechanistic  pyrolysis 
model  for  cellulose  pyrolysis  provided  rate  equations  for  volatiles  and  char  which  reduced  to  a 
single  first-order  rate  law  for  the  isothermal  mass  by  making  the  reasonable  assumptions:  (1) 
primary-bond  dissociation  is  the  rate-limiting  first  step,  (2)  mass  loss  proceeds  through  an  active 
intermediate  which  is  a  stationary  state,  (3)  primary  gas  and  char  production  are  competing 
processes  whose  magnitude  and  rate  constants  are  large  compared  to  successive  products  and  rate 
processes,  and  (4)  conditions  are  anaerobic.  Exact  results  were  obtained  for  the  isothermal  mass 
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loss  history  of  a  char  forming  polymer  which  included  an  equilibrium  char  yield  whose  value 
depends  only  on  temperature  through  the  relative  rates  of  gas  and  char  formation.  In  the  present 
analysis,  the  crossover  temperature  at  which  rates  of  gas  and  char  formation  are  equal  corresponds 
to  the  temperature  at  which  the  equilibrium  char  yield  is  50%. 

Comparison  of  the  model  prediction  with  the  measured  mass  loss  history  of  a  phenolic  triazine 
thermosetting  resin  in  temperature  scanning  thermogravimetric  experiments  showed  excellent 
agreement  for  the  first  stage  of  thermal  decomposition  to  volatiles  and  char.  An  approximate 
expression  for  the  maximum  mass  loss  rate  at  constant  heating  rate  was  derived  which  is  in 
qualitative  agreement  with  thermogravimetric  data  for  charring  and  noncharring  polymers.  The 
simple  analytic  solution  for  peak  mass  loss  rate  during  steady  heating  shows  that  the  maximum  fuel 
generation  rate  decreases  linearly  with  char  yield  and  is  inversely  proportional  to  the  square  of  the 
(rate  dependent)  peak  decomposition  temperature.  Thus,  the  analytic  result  for  mass  loss  in  a 
relatively  complex  polymer  thermal  decomposition  process  is  physically  based  with  temperature 
and  heating  rate  dependent  product  yields  which  are  important  in  fire  modeling  (e.g.,  volatile  and 
char  formation).  Although  outside  the  scope  of  the  present  work,  the  incorporation  of 
thermoxidative  processes  in  the  global  thermal  degradation  scheme  (equations  6-8)  allows  for  a 
more  detailed  analysis  of  the  mass  loss  rate  of  polymers  heated  under  aerobic  conditions. 
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